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Abstract 

A comparative study of the numerical renormalization group and non- 
crossing approximation results for the spectral functions of the U = oo An- 
derson impurity model is carried out. The non-crossing approximation is the 
simplest conserving approximation and has led to useful insights into strongly 
correlated models of magnetic impurities. At low energies and temperatures 
the method is known to be inaccurate for dynamical properties due to the 
appearance of singularities in the physical Green's functions. The problems 
in developing alternative reliable theories for dynamical properties have made 
it difficult to quantify these inaccuracies. In this paper we show, by direct 
comparison with essentially exact numerical renormalization group calcula- 
tions for the auxiliary and physical particle spectral functions, that the main 
source of error in the non-crossing approximation is in the lack of vertex cor- 
rections in the convolution formulae for physical Green's functions. We show 
that the dynamics of the auxiliary particles within NCA is essentially cor- 
rect for a large parameter region, including the physically interesting Kondo 
regime, for all energy scales down to To, the low energy scale of the model 
and often well below this scale. Despite the satisfactory description of the 



auxiliary particle dynamics, the physical spectral functions are not obtained 

accurately on scales ~ To. Our results suggest that self-consistent conserving 

approximations which include vertex terms may provide a highly accurate 

way of dealing with strongly correlated systems at low temperatures. 
PACS numbers: 71.27.+a,71.28.+d,72.20.Hr 
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I. INTRODUCTION 



In recent years the physics of strongly correlated Fermi systems has attracted wide in- 
terest, in particular in the context of high T c superconductors Q and their normal state 
properties 0. Despite intense efforts, systematic and controlled theoretical methods for 
dealing with such systems are still lacking. At the heart of the problem is the effect of a 
strong on-site Coulomb repulsion U on fermions living on a lattice, in particular in d = 2 
dimensions. At low energy such models can be mapped onto effective Hamiltonians with the 
constraint of no double occupancy of any lattice sites. The projection onto the correspond- 
ing subspace of Hilbert space is presumably responsible for a number of unusual properties, 
most remarkably the spin-charge separation and non-Fermi liquid behaviour. Unfortunately 
the local constraints on the site occupancy are difficult to handle. One possible formulation 
uses auxiliary particles such as slave-bosons (6) and pseudo-fermions (f a ) to represent an 
electron as a composite particle consisting of a pseudo-fermion particle and a slave-boson 
hole {fib). The local constraint is holonomic in this formulation, and is given by the condi- 
tion that the sum of the occupation numbers of pseudo-fermions and slave-bosons is equal 
to unity at any lattice site, J2a n f(r + n b = 1. The constraint is closely related to the local 
gauge symmetry of the system with respect to simultaneous U(l) gauge transformations of 
b and f a . It may be shown that starting with certain mean field theories the constraint and 
its consequences lead to the appearance of longitudinal and transverse fictitious gauge fields 
coupling equally to pseudo-fermions and slave-bosons [|J. The gauge fields restore the local 
symmetry broken in mean field theory. 

However, it is not clear that the slave-boson mean field theories are a good starting point 
(except for the cases where a physical phase transition into a magnetic or superconducting 
phase takes place). Alternatively one may guarantee the local gauge invariance of approx- 
imations by deriving those from a generating functional ||. The local gauge invariance 
ensures conservation of the local occupation at each lattice site (the "local charge") in time. 
In order to effect the actual projection, it is still necessary to fix the occupation numbers 
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at a given time locally. So far this has only been achieved for impurity models, e.g. the 
Anderson model of a magnetic impurity in a metallic host 

Even if the projection is done exactly, an equally important question is the selection of 
the dominating contributions in perturbation theory in the hybridization (or in the hopping 
integrals and exchange interaction in the case of lattice models). In this paper we address 
this question for the infinite U Anderson impurity model (degeneracy N = 2). We employ 
the numerical renormalization group (NRG) to calculate the slave-boson and pseudo-fermion 
spectral functions and compare the results with those of the simplest fully projected con- 
serving approximation, the so called "non-crossing approximation" (NCA) @JT^] • While the 
NCA is known to give excellent results for large degeneracy N and not too low temperatures, 
including the cross-over from low (T « T K ) to high temperature (T >> T K , where T K is 
the Kondo temperature) 0,0, pronounced deviations from exactly known results appear 
at iV = 2 |ynj]. The comparisons will allow us to pinpoint the deficiencies of the NCA and 
to identify possible improvements. 

In earlier publications ]T2"| , P~5[ we have presented results of an NRG calculation of the 
slave-boson and pseudo-fermion spectral functions at zero temperature. We found these 
functions to be infrared divergent at threshold, with critical exponents dependent on the 
d-level occupancy and given by simple expressions identical to the well known X-ray 
absorption threshold exponents. 

The paper is organized as follows: In Section II we formulate the model, describe how 
we implement the numerical renormalization group for studying auxiliary spectral functions 
of this model and describe the NCA within the framework of conserving approximations. 
Section III describes our results for spectral functions calculated with the above two meth- 
ods. In Section IV we summarize our main results. Some details of the auxiliary particle 
technique, which make it suitable for an effective evaluation at the lowest temperatures, are 
discussed in the appendices. 
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II. FORMULATION 



A. Model and Application of the Renormalization Group 

The Anderson model of an impurity d electron state hybridizing with the conduction band 
with infinitely strong Coulomb repulsion in the d-level in auxiliary particle representation 
takes the form 

H = H c + e d J2flf* + Vj2(^fa + h.c), (1) 

(7 (7 

where H c = Yl,ko e kc\ a Cko is the conduction electron kinetic energy and c 0cr = J2k c ka an- 
nihilates a conduction electron with spin a at the impurity site 0. The Hilbert space 
consists of disjoint subspaces characterized by the conserved auxiliary particle number 
Q = tfb + Yla fa-fen Q — 0,1,2,.... The physical subspace is defined by the constraint 



Q = 1. Following Wilson |14| we (i) linearize the spectrum of H c about the Fermi energy 
6fc — > k, (ii) introduce a logarithmic mesh of k points k n = A~ n and (iii) perform a unitary 
transformation of the Ck a such that Co a is the first operator in the new basis and H c takes 
the form of a tight-binding Hamiltonian in k-space, 

oo 

H c = E E£nA- n/2 ti + i<x<W + h.c), (2) 

71=0 V 

with £ n — > (1 + A _1 )/2 for n >> 1. These steps are explained in detail in |14| and can be 
taken over for the present model without change. 

The Hamiltonian ([!]) together with the discretized form of the kinetic energy (0) in the 
new basis is now diagonalized by the following iterative process: (i) one defines a sequence 
of finite size Hamiltonians H N by replacing H c in (g) by H C N = Y,n=o E CT Cn^' n/2 (cl +la c na + 
h.c); (ii) starting from H = eaJ2a /]/(T + ^E ff ( c L^/ff + ^ c -)i each successive hopping may 
be considered as a perturbation on the previous Hamiltonian; (iii) the Hamiltonians are 
scaled such that the energy spacing remains the same. This defines a renormalization group 
transformation H N+1 = k l / 2 H N + ^nUA^n+i^No- + h.c.) — Eq,n+i, with Eq,n+i chosen 
so that the ground state energy of Hn+i is zero. The Hamiltonians Hn are diagonalized 



within subspaces of well defined Q, N e , S, S z where N e is the total fermion number, S the 
total spin and S z the z-component of total spin. The use of these conserved quantities leads 
to significant reductions in the size of matrices to be diagonalized, however the dimension 
of Hn grows as 4^ and it is necessary to truncate the higher energy states for N > 7. 
Approximately 1 /4 of the states generated at each iteration are retained in the calculations 
and this constitutes a surprisingly accurate approximation for the present model. This 
accuracy is evidenced by the fact that various exact relations, such as the Friedel sum rule, 
which relates the impurity spectral density at the Fermi level to the local level occupancy, 
are satisfied to a high degree of accuracy [|15|]. In addition, the number of states retained 



per iteration, N st , which is a free parameter, can be varied from 250 to 2000 states without 
any significant change in the results presented here. This, together with the almost cut-off 
independent results for A > 1.5 indicates the convergence and accuracy of the method. 

B. NRG Calculation of Auxiliary Spectral Functions 

Within the framework of the NRG it is natural to represent the Green functions in an 
eigenbasis of the Hamiltonian. In the enlarged Hilbert space the retarded propagators for 
pseudo-fermions and slave-bosons are defined by 

G fa (co + zO, T,\) = ^~ E \Mlj 2 (e-^ +XQ ^ + e~^ +x ( Q+1 ») / (to + z'O - A — (e m - e n )), 

Q,m,n 

G b (oo + iO, T, A) = -J— X \M b mn \ 2 (e-^ +x ® - e -^ +x ^)/(u + iO - X - (e m - e n )), 



Z GC q 



where, for each Q, e m , e n are eigenvalues of the subspaces Q + 1, Q, where Z GC (T, A) is the 
grand-canonical partition function, 

Z GC (T, A) = X; e-^ +AQ ) = Z Q=0 + e-? x Z Q=1 + ... (3) 

Q,m 

and M% n =< Q + l,m\0^\Q, n > with O = f a ,b are the many-body matrix elements 
for pseudo-fermions and slave-bosons, respectively. The constrained propagators are then 
obtained as 
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G ftb (u + zO, T, A - oo) = E \Mt h n \ 2 e-^/{uj + zO - (e m - e n )), (4) 

Z Q=0 m ,n 

where the frequency has been shifted as w — > uo + A. We are interested in 
the following projected spectral functions (see Appendix C for details)) A^ b (u,T) = 
-lim^oopm G f>b (uj,T,X)/ir) and Aj jb (u,T) = -lim^oo -§^[lm G f:b (uj,T,X)/n], 



A+^iu, T) = 1 £ | < 1, m \0 ] |0, n > | V^°.«% - (e 1>m - e ,„)), (5) 



m,n 



Aj >b (u;,T) = — £ | < l,m\0%n> \ 2 e~^6(u - (e ljm - e 0l „)), (6) 

In the definition of Aj b Zc = Zq = i is introduced in order to obtain a well defined zero 
temperature limitQ. At zero temperature the spectral functions reduce to 

A+ b (u, T = 0) = £ I < 1, m\O^ > Mu: + E% s =0 - e 1>m ), (7) 



z Q=0 (o) 
1 

Zcifi) Z Q=0 (0) 

Here |$o > is the groundstate of the Q = subspace of non-interacting conduction electrons 



Aj >b (u;, T = 0) = - ^ E I < $ i|O f |0, n > + e 0) „ - JSgSj. (8) 



and |l,m > are the excited states of the Q — 1 subspace of the interacting system, -Eq£ 
and ei )m are the corresponding energy eigenvalues. The spectral functions A^ b {T = 0) 
vanish identically below the threshold E = EqZ^ — Eq= . Similarly, in the expression 
for Aj b (T = 0), |$i > is the groundstate of the interacting system (Q — 1 subspace) 
and \0,n > are the excited states of the non-interacting conduction electron system, Eq^ 
and eo, n the corresponding energy eigenvalues. These spectral functions vanish above the 
threshold energy Eq. From ( |B4"1 - [B5|) and ([?]) we see that the T = spectral functions A^ b 



satisfy the sum rule 

r+oo 

A+ b {uj,0)duj=l (9) 



1 The functions , A b defined here correspond to the A, B functions, respectively in Ref. 0, 
and similarly the functions Aj , A b correspond to the a and b functions. 
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In practice, within the NRG technique, it is the excitation energies from the respective 
Q = or Q = 1 groundstates which are calculated. Hence, it is convenient to set the 
threshold energy, Eq, to zero. The T = projected auxiliary spectral functions then have 
divergences at zero energy. As described in Appendix D it is also possible to formulate 
the NCA equations of section C such that the NCA spectral functions exhibit their T = 
divergences at zero energy. This makes it easier to compare with the NRG results and in ad- 
dition improves the accuracy of the NCA solution at the lowest temperatures. The reference 
energy scale for the physical d electron Green's function, Gda(^,T) =« b'f a ; bf\ »q=i, 

G^uj + *0, T) = ]T \M d m J\e~^ + e-*»)/(o; + *0 - (e m - e n )), (10) 

C m,n 

is independent of E and is set by the Fermi level ep = 0. 

The matrix elements M^ b n of the pseudo-particle operators /t, in (||-§D are calculated 
recursively using the formulae given in Appendix A. Similar formulae apply to the matrix 
elements Mf nn for the physical d electron Green's function. For each iteration step, they 
are substituted together with the energy eigenvalues into to give the T = spectral 

functions A N j b (u>). In principle if all states up to stage N were retained, Hn would describe 
excitations on all energy scales from the band edge D = 1 down to the lowest energy scale 
present in H^, i.e. u>n- Due to the elimination of higher energy states at each step, the 
actual range of excitations in is restricted to uj^ < uj < Ku^, where K « 7 for A w 2 
retaining 500 states per iteration. Thus at step N, the spectral functions are calculated at an 
excitation energy uj w 2un in the above range. The delta functions in (|7|-§) are broadened 
with Gaussians of width ~ appropriate to the energy level structure of H N [|15| . 



Finally we note that the constraint of no double occupancy, Q = 1, is implemented 
exactly within the NRG calculations since Q is a conserved quantum number. Appendix C 
describes in detail the implementation of this constraint for the conserving approximation 
of following section. 
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C. Conserving approximation: NCA 



For small hybridization V compared to the impurity level e^, a perturbation expansion in 
terms of V seems reasonable. It is useful to express the auxiliary particle Green's functions 
for given "chemical potential" —A defined in ( |CT| ) in terms of self-energies E&(u;,T, A) and 
Tif(u>, T, A) as 

G fa (u, T, A) = [u - A - e d - E f (u, T, A)]" 1 (11) 
G b (u, T,X) = [u-X- E b (oj, T, A)]- 1 (12) 

where u takes the values of fermionic or bosonic Matsubara frequencies. The projected 
spectral functions are obtained by shifting the frequency uo — > uo + A and taking the limit 
A — > oo, as described in Appendix C. The local Green's function of conduction electrons 
G C(7 (uj, T, A) is given likewise in terms of the self energy X ccr (c<j, T, A) as 

G ca (uj,T,X) = [(G^y 1 - S CCT (u;, T, A)] -1 (13) 

where G C(T {uj) = J deN(e)/(u — e) is the local Green's function of the bare conduction band, 
with N(e) the density of states. There is an exact relation between the d-electron Green's 
function and E c : 

G " ( "' T ' A) = ^i-Zt^c^r (14) 

Thus, the principal problem remaining in the evaluation of the theory is the proper choice 
of approximation in which the three self-energies are calculated. The slave-boson Hamilto- 
nian ([!]) is invariant under two independent U(l) gauge transformations: (1) simultaneous 
transformation of the pseudo-fermion and the slave-boson operators and (2) simultaneous 
transformation of the pseudo-fermion and the conduction electron operators. These symme- 
tries correspond to the conservation of the auxiliary particle number and the total number 
of fermions, respectively: 

-£fj a + tfb = const. (15) 
H flU + X c L°k* = const. (16) 

a ka 
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The Green's functions Gf, Gb are not invariant under local (in time) gauge transformations, 
but we must require Gd to be invariant. Constructing a gauge invariant approximation is, 
a priori, a non-trivial task, since the d electron number rid is not a conserved quantity. 
Yet, it may be shown that a gauge invariant approximation obeying both conservation laws 
(|T5D, QTjSp is constructed by deriving all three self-energies from one generating functional 
$ : Eyv = 5§/8Gf a , S fe = 5$ /8G b , S C(T = 5§/5G ca . Eq. fllj) then provides the rule for a 
gauge invariant approximation for Gd- The functional <3> is given in terms of closed skeleton 
diagrams with suitable combinatorial factors. The lowest order contribution to $ is of second 
order in V, $ = — V 2 ^ J2u S e G f a (u)G ca (e)Gb(uj — e) (Fig. 1 a). Functional differentiation 
yields the self-energies 

S /<T ( W ) = -V 2 ± £ CUe)G 6 (w - 6) (17) 
£ 6 (fi) = V 2 i £ G /(7 (ft + e ) G -( e ) ( 18 ) 

P e 

= 2 G fff (e)G b (e - u). (19) 

P e 

After the transformation uj — > u + \, A — > oo, one finds explicitly 

£ /(T (u;,T) = V 2 ]T(1 - f(e k ))G b (u - e k ) (20) 
£ 6 (u;, T) = V 2 J2 f(^)G fa (u + e k ), (21) 

k.a 

E„(u, T, A -> oo) = ^ 2 e^ A | dee-^[G /CT (e + w)A+(e) - A+(e)G b (e - w)] (22) 

and G/^u;) -1 = cu— e^— £/o-(co>), G^la;) -1 = u—T,b(uj), where /(e^) is the Fermi function. The 
physical d-electron Green's function is obtained from the limiting procedure Gd a (yi,T) = 
lim^-,00 e /3A G ( fo-(a;, T, A) (compare Appendix C). As E co .(u;, T, A) ~ e~ l3X for A — * oo, the 
self-energy corrections to G c and to the denominator of Gd vanish, and 

GW(^,T) = i- lim e^E^T.A). (23) 

The impurity spectral density follows as, 
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^ CA (u, T) = --Im G da (uj + zO, T) 

7T 
r+oo 



/+oo 
de[A+{e + w, T)A 6 -(e, T) + Aj(e, T)A+(e - w, T)], (24) 
-oo 



which at T = reduces to, 



^, T = 0) = 6(cj) / deA+(e + u, T = 0)A^(e, T = 0) 

+ e(-w) / rfeA7(e,T = 0)^(e-c;,T = 0). (25) 

JEq+oj 

The above approximation is known as the "non-crossing approximation" (NCA) because it 
does not include any diagrammatical contributions with crossed conduction electron lines 

In general the impurity electron Green's function may be expressed with the help of a 
vertex function A(e, u) as (Fig. 1 b) 

G da (cu) = - lim e^J2G f a(e)G b {e-u)A(e,Lo) (26) 

e 

As shown below, there is reason to expect that the vertex function A plays an important 
role. 



III. RESULTS 

The NRG calculations were performed for A = 2, keeping 250 states per iteration for each 
subspace (Q = 0, 1). The hybridization strength A = it V 2 p(ep) = vr V 2 /2D was chosen to 
be 0.01-D with the half-bandwidth D — 1. Several values of the local level position e d were 
chosen in order to characterize the behaviour of the spectral densities in the various regimes. 
The NRG spectral functions were evaluated at T = and the NCA spectral functions were 
evaluated following appendix D for the same set of parameters and for temperatures down 
to T = 10~ 6 D << To, where To is the low energy scale of the model. This was sufficiently 
low to allow comparison with the T = NRG results over most of the interesting energy 
range. We define Tq to be the Kondo temperature, 



11 



k B T K = DJ-e~^ 2A (27) 

in the Kondo regime e^/A < —2, A in the mixed valent regime |erf/A| < 1 and in the 
empty orbital regime e^/A >> 1. 

A. Threshold behaviour of the NRG auxiliary spectral functions 

The T = auxiliary spectral functions diverge at the threshold Eq as shown in Fig. (0- 
^|). This behaviour may be understood as a result of the orthogonality catastrophe theorem 
. To see this more clearly we re-formulate the spectral densities in in the following 

way, 

A} tb {u, T = 0) = — !— £ | < 1, m|$ > | 2 5(^ - ei,m), (28) 

Z Q=0 m 

A7 fe (cu, T = 0) = — !— £ I < <£ x |0, n > | 2 5(cu + eo, n ). (29) 

Z Q=0 n 

In the above, |<3>o >= Ot|$ >> with O = f a ,b, represents the non-interacting (U = 0) 
groundstate with 1(0 = fl) or (O = tf) local electrons present. Similarly, |$i >= 0|$i >, 
with O = / CT , b, represents the interacting (U = oo) groundstate with 1 (O = /j) or 
(0 = 6^) local electrons present. In this formulation, we see that A^ b measures the overlap 
density between the groundstate of the non-interacting (U = 0) band electrons with or 1 
local electrons present with the excited states of the interacting U = oo Hamiltonian. Sim- 
ilarly, Aj b measures the overlap density between the groundstate of the interacting U = oo 
Hamiltonian with or 1 local electrons present and the excited states of the non-interacting 
band electrons. This interpretation is identical to that for the core level spectral functions 
in the X-ray problem. The analogy is useful but requires care since the matrix elements 
in (0-§D are no longer between two non-interacting systems as in the X-ray problem. This 
leads in particular to a new energy scale, To, for the onset of the asymptotic power law 
behaviour, which is TV, A or in the Kondo, mixed valent and empty orbital regimes, 
respectively. We find that it is only in the Fermi liquid regime, \uo — E \ << T , that the 
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power law behaviour is well characterized. The approach to this asymptotic power law is 
faster for the boson spectral functions than for the fermion spectral functions in all cases. 
We note that within the NCA, the approach to the threshold behaviour with n d independent 
exponents a/ = 1/3 and a& = 2/3 is also only asymptotic and requires in particular going 
down to temperatures T < 1CT 2 T in order to see these exponents. 

The threshold exponents for the slave-boson and pseudo-fermion spectral functions were 
extracted by numerically differentiating the spectral functions. Typically, well defined ex- 
ponents can be extracted only for energy scales \u — Eq\ < lCT 2 To. The exponents are 
shown in Fig. (|5]) and Table | as a function of n<2, the local level occupancy at T = 0. The 
latter, shown in Fig. was calculated by evaluating n^(T) from the partition function at 
a sequence of decreasing temperatures T/v ~ A"^ -1 ^ 2 and then taking the limit T — > 0. 

Remarkably, the threshold exponents turn out to be the usual photoemission and ab- 
sorption exponents for the X-ray problem and are given in terms of the conduction electron 
phase-shift at the Fermi level [EJ, 5 a = ^(e^), by 

«/ = - d -f = 2^-E(-) 2 (30) 

2 (7 

2 „ n 



where the last equations on the RHS of (pO^pT]) follow from the Friedel sum rule, 5 a = it n d /2. 
These results are clearly illustrated in Fig. (|5|) where the functions — n 2 d /2 and 1 — n 2 d /2 
are plotted against n d together with the exponents a deduced from the spectral functions. 
The exponents acf^ agree with the RHS of fl30|-|3~T|) to 3 significant figures in nearly all cases 



and are the same below and above the threshold, 

A% h = a% b \u - E \-^» (32) 
A qualitative argument based on charge neutrality considerations has been given for the 



above form of the exponents \12\. We note that the same functional form of the exponents on 



the phase shift (|3T)|-pT|) is also found in the spinless model with constraint jTT3| in agreement 



with exact analytic results JT^|. The single phase shift in this case is given by 6 = nn^ 
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. An rid dependent exponent was also found by Read in considering how Gaussian and 
higher order corrections restore the gauge symmetry broken by the slave-boson mean field 
theory [|Tj| . Generalizing the above threshold exponents to the N-fold degenerate model 
we have a f = 25 m /'n - Eml'W 71 ") 2 = 2n d/N - n 2 d /N and a b = 1 - E m (^mA) 2 = 1 - n 2 d /N 
where m labels the scattering channels. The same exponents were found for the N— fold 
degenerate Anderson impurity model in perturbative calculations to order V A adding to 
the plausibility of the above generalization. The above conjecture for af^ is in disagreement 
with recent results obtained in the limit rid —* 1 and in the large N expansion including order 

^ 1 AT 2 

1/N 2 |p0| , a/ = jyTT~r ) O-b = jvqr—r • These results were obtained using a perturbative 

iV 2 " N 2 

renormalization group technique, which we do not expect to be as accurate as the non- 
perturbative numerical scheme used here. We see from our results, generalized to arbitrary 
N, that the NCA exponents 1/(JV+1) = 1/N+0(1/N 2 ) andiV/(iV+l) = 1-1/N+0(1/N 2 ) 
are correct only in the limit rid — > 1 and N — > oo (or in the trivial limit rid 0). Away 
from this limit, vertex corrections in the auxiliary Green's functions, absent in the NCA, 
are therefore important in determining the correct threshold exponents. The expressions for 
the threshold exponents of auxiliary particle propagators in terms of X-ray photoemission 
exponents appears to be a general property of several impurity models exhibiting Fermi 
liquid fixed points. 



B. NRG auxiliary spectral functions at higher energies 

At higher frequencies the following features are observed in the T = spectral functions. 
In the Kondo regime, Fig. fl2|), there is a peak in the slave-boson spectral function A£ at 
uj = | e<2 1 and a much less pronounced feature in the corresponding pseudo-fermion spectral 
function Af. As is raised through the Fermi level from below the peak in A^ at u — | 
becomes less pronounced and almost disappears in the mixed valent e^/A ~ 0, Fig. 
and empty orbital regimes e^/A > 1 Fig. (|4]). In addition, its position is renormalized 
above the bare value \ed\- At the same time the small feature at positive energies in the 
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pseudo-fermion spectral function Aj develops into a well defined peak in the mixed valent 
and empty orbital regimes. The "— " spectral functions Aj b (u, T = 0) exhibit monotonic 
behaviour for all parameter regimes. 

C. Comparison of NRG and NCA auxiliary spectral functions 

In comparing NRG and NCA spectral functions, three energy regimes should be distin- 
guished: 

(I) asymptotically low energy regime, \u — E Q \/T « 1, 

(II) crossover regime, \u — E \/T ~ 1, 

(III) high energy regime, \lj — E \/T Q » 1. 

The energy range (I) corresponds to that discussed in the section on threshold exponents. 
Here we discuss the energy ranges between (I) and (II), and between (II) and (III). In Fig. 
(0H1) the same qualitative trends described in the previous section for the NRG auxiliary 
spectral functions can be seen in the corresponding NCA solutions. The NCA results for 
a finite but very low temperature T = 10~ 6 D are compared to the corresponding T = 
NRG results. Surprisingly good quantitative agreement is seen in the slave-boson spectral 
function A£ above the threshold down to energy scales well below T in all regimes. The 
agreement is particularly good in the Kondo regime for —4 < e^/A < —2, where it extends 
down to 10~ 2 T (e.g., Fig. |7|a). The spectral function below the threshold also shows 
good agreement with the NRG result in the Kondo regime with decreasing agreement in 
the mixed valent e^/A ~ and empty orbital regimes e^/A >> 1. Turning now to the 
pseudo-fermion spectral functions we see that there is again good quantitative agreement 
between NCA and NRG for above the threshold and for all energy scales down to 
To. This is true in all regimes. Below the threshold the agreement for the A J spectral 
function even extends to well below To, except in the Kondo regime for e^/A < —4 where 
we could not obtain quantitative agreement except in the region 10 _1 < ui < 10 1 . From 

15 



these comparisons we see that the most serious differences, as far as low energy behaviour 



is concerned, between the NCA and NRG auxiliary spectral functions are in At for u < T 



in the Kondo case and in in the mixed valent and empty orbital cases. The latter we 
attribute to the inaccuracy of the NCA in the energy range between (I) and (II) in the 
mixed valent and empty orbital regimes to be described further in the section on impurity 
spectral densities. A more interesting discrepancy which arises from these comparisons is 
the former. The functions A^ b in NCA are related by self-consistency equations (|20| |2T|) 



derived by second order perturbation theory in V. As pointed out in p2| , coherent spin 
flip processes considered to be responsible for the Kondo resonance are not included in the 
NCA. There are reasons to expect that the self-consistent T-matrix approximation proposed 



m 



22 1 will capture the essential contributions. Recent calculations within the conserving 



self-consistent T-matrix approximation indicate that such improvements do indeed arise 



A different extension of NCA, called "post-NCA" has been proposed recently on the 
basis of a expansion of the N-orbital model |2T|. It represents a self-consistent scheme 
including vertex renormalizations, which is exact to order 1/N 2 . 

D. Comparison of NRG and NCA impurity spectral functions 

In the previous section we noted that the NCA auxiliary spectral functions were surpris- 
ingly close to the NRG ones for energies down to at least u = Tq and typically they were 
even quantitatively accurate down to uj « Tq. Improvements are primarily important in 
the auxiliary spectral functions in two areas, to restore the correct behaviour of the pseudo- 
fermion spectral function below To and to recover the exact threshold exponents given 
by the NRG. We now turn to the comparisons for the impurity spectral function and discuss 
the role of vertex corrections on the dynamics of the physical electrons. The impurity spec- 
tral densities are shown in Fig. (|IO|-|Il2]) where in addition to the NCA, p^M,and NRG 
spectral functions, p^ RG (tu), we also show the impurity spectral function, p^ RG ~ , obtained 
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by convoluting the NRG auxiliary spectral functions as in (|2q) but without the vertex part, 
i.e. 

/+oo 
de[A+(e + lu, T)A»(e, T) + Aj(e, T)A+(e - w, T)). (33) 
-oo 

In the Kondo regime (Fig. 110 a-c), the impurity spectral density shows two peaks, a charge 
fluctuation peak close to the local level position e^, and a many-body Kondo resonance 
at the Fermi level. The charge fluctuation peak is broader in the NRG case, a result of 
using a logarithmic discretization for the conduction band which leads to lower resolution 
at higher energies. This is not a fundamental problem with the NRG, and the resolution of 
the method at higher energies can be improved by reducing the discretization parameter A 
for the relevant iterations covering the high energy scales. We have explicitly checked that 
the width and height of the many-body resonance at the Fermi level, where the NRG gives 
the highest resolution, are unaffected by the broadening used to smooth the delta functions 
in the discrete spectra (unless the broadening is made too small which will result in uneven 
spectra). From Fig. pH|b-c we see that the NCA gives a Kondo resonance which is too broad 
and too high on the BIS (u > 0) side for energies below ~ 5T . On the PES side (u < 0), 
i.e. for {jj < 0, the agreement with the NRG is better. 

The disagreement between the NCA and NRG impurity spectral function for energies 
below 5T is due primarily to the absence of the vertex part in the convolution formula for 
the NCA impurity spectral density (0). This is seen from the good agreement between 
Pd RG ~ an d Pd° A m the range T < uj < 5T which indicates that the NCA auxiliary 
spectral functions are sufficiently accurate in this range and that therefore the difference 
between p^ RG and p^ CA must be mainly due to the neglect of the vertex part in ^j). Similar 
comparisons for the other Kondo cases e^/A = —5,-3 support the same conclusion. In 
Table | we also list p^ RG {uj = 0) and p^ CA (cu = 0, T = 1CT 6 -D) together with the respective 
relative deviations from the exact Friedel sum rule result, Pd a (uj = 0) = sin 2 (7rn ( i/2)/7rA. 
The NCA result in the Kondo regime appears reasonable because the singular behaviour 



of the impurity spectral density at T = JTTJ] is removed by our small finite temperature 
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T = 10 6 D. The exact NCA result for p^ CA {oj = 0, T = 0) is actually much worse [11 . 



We also show results for the mixed valent, Fig. QTT|), and empty orbital regimes, Fig. 
(|T2|). There is good agreement between p^ RG , Pd RG ~ an d Pd° A f° r the high energy parts of 
the renormalized resonant level uj > T , but the incorrect low energy behaviour of the NCA 
result for in the above regimes and the neglect of the vertex part in ( p^) and (|33| ) makes 
both p^*"" 4 and p^ RG ~ deviate from the exact NRG result at low energies. The resonant 
level is approximately a Lorentzian of width A, and the small asymmetric broadening in 
the NRG curves is due to the logarithmic discretization. An improved description of high 
energies could be obtained in both NRG and NCA, if required. Within NRG it is possible 
to focus on high energies explicitly by using a smaller discretization parameter, A, for the 
first few iterations. Within NCA higher energies are easily resolved by using a finer grid to 
solve the integral equations at these energies. 

Having discussed the effect of vertex corrections on the impurity spectral function at 
intermediate u ~ T and higher energies u > T in the different parameter regimes, we now 
discuss the limit uj «Tq. In this limit, evaluating the impurity spectral density, without 
vertex corrections, i.e. for A = 1, using (|30|-p^) in fl26|) gives 



p da (uo -> 0+) = a+a^uj^ a f- ab B(l - a f , 1 + a b ) ~ \u\- nd{l - nd) . (34) 

where B is the Beta function. The exact result at to = is given by the Friedel sum rule, 
Pda{u = 0) = sin 2 (7mrf/2)/7rA, so we conclude that the vertex corrections neglected in (|3~4"D 
are singular at low energies, i.e. close to the threshold, and lead to a singularity in pd a at the 
Fermi level which cancels that in (0). Similar vertex corrections appear in the calculation 
of other physical quantities such as the dynamic spin susceptibility. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we have made a comparison of the spectral functions of the U — oo Anderson 
model as calculated within the simplest self-consistent conserving approximation, the NCA, 
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and within the NRG. At high energies u » To we found good quantitative agreement for the 
auxiliary and physical spectral functions calculated within the two methods in all parameter 
regimes. Some small discrepancies in the shape of high energy peaks could be attributed 
to the logarithmic discretization used in the NRG which tends to give lower resolution and 
slight asymmetries to high energy peaks. At lower energies we found good quantitative 
agreement for both slave-boson and pseudo-fermion spectral functions down to at least 
To. In the Kondo regime the agreement between the NCA and NRG slave-boson spectral 
function extended to well below the low energy scale To = Tk- Despite the accuracy of 
the NCA auxiliary spectral functions down to Tk we noted that the impurity spectral density 
deviated from the essentially exact NRG result on energy scales up to 5T K . The source of 
this discrepancy was traced directly to the lack of the vertex part in the NCA expression 
for the impurity spectral density, thus showing that vertex corrections are required for the 
physical Green's functions even when the auxiliary particle dynamics appears to be correctly 
described within the NCA (for energies down to T ). In the Fermi liquid regime u << T , 
the NCA gives results for the impurity spectral density in disagreement with exact results 
from Fermi liquid theory. It also gives the incorrect threshold exponents aif^ for the auxiliary 
particle Green's functions. We emphasize, however, that although vertex corrections restore 
the Fermi liquid behaviour in the physical spectral functions and the correct low energy 
asymptotic behaviour of the auxiliary particle Green's functions, our results show that they 
are also needed for a correct description of the physical Green's functions at much higher 
energies u ~ T . This is the main conclusion of the comparisons we made. 

We are in the process of including vertex corrections within a self-consistent theory |22| . 
thus going beyond the lowest order conserving approximation. The availability of accurate 
results for dynamic properties via the NRG for impurity models makes the latter a natural 
testing ground for developing such approximation schemes. These schemes could then be 
extended to study lattice models of strongly correlated electrons, for which, at present, 
numerical renormalization group methods are not as well developed as for impurity models. 
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APPENDIX A: MATRIX ELEMENTS FOR AUXILIARY PARTICLE 

OPERATORS 

In this appendix we give the expressions for the matrix elements of the auxiliary 
particle operators, required for the calculation of the spectral functions. The states, 
\Q, N e , S, S z ,r >at, of the Hamiltonian Hn (which includes the orbitals f a , co a , . . . , cno) 
are labeled by the quantum numbers Q, N e , S, S z and an index r, where Q = J2a ft fa + b^b 
is the number of auxiliary particles, N e is the total number of electrons, S the total spin 
and S z the z-component of the total spin and the index r distinguishes states with the 
same conserved quantum numbers. A product basis set \Q, N e , S, S z , ri >n for the subspace 
(Q, N e , S, S z ) of H N in terms of eigenstates of H N _i and states |0 >, | |>, | |>, | f|> of the 
orbital cjy a is defined by 

\Q, N e , S, S z , r,l> N = \ Q N e ,S, S z , r > N ^ |0 >, (Al) 
\Q, N e , S, S z , r, 2 > N = fi^K\Q, N e -l,S-±,S g -±r >at_i | T> 

+ \[^^\Qi N e -l,S-±S, + ±r >at_i | |>, (A2) 
\Q,N e ,S,S z ,r,3> N =- ] j S ~ S S ^ f |g, N e -l,S+±,S g -±,r >JV -i I T> 

+ \j S 2 S S l2 1 |g ' Ne ~ 1j S + h Sz + \ >N ~ X 1 i>5 (A3) 
\Q, N e , S, S z , r, 4 >jv = \Q, N e - 2, S, S z , r > N ^ \ ||> (A4) 

The reduced matrix elements n < Q, N e , S, r\ \f^\ \Q — l,A^ e — 1,S ± 1/2, s > N and jv < 
Q, N e , S, r\ \Q — 1, A^g, S, s >n are calculated recursively following the recursive calculation 
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of the matrix elements n < Q, N e , S, r||cjy||Q, N e — 1, S ± 1/2, s >at required for setting up 
the Hamiltonian H^. Details of the latter can be found in ||14)| . We follow the notation 



of 0] and denote by Uq NeS (ri,p),p = 1, . . . , Rq NeS the matrix of eigenvectors of the 
subspace (Q, N e S, S z ) of Hn where Rq^s * s the dimensionality of this subspace and r and 
z = 1, . . . , 4 label the product state basis which is related to the diagonal basis by the unitary 
transformation 

\Q,N e ,S,S z ,p >N=J2 U Q,Ne,s(p,ri)\r > N ^ \i > (A5) 

r,i 

with |r >jv-i \i > denoting one of the four product states (|7n|- |A4| ) defined above. Defining, 



K% N Ar, r') = N < Q, N e , S, r\\p\\Q - 1, N e - 1, S ± i r' >7V (A6) 



MSI Ar, r') = N < Q, N e , S,r\\tf\\Q ~ 1, N e , S, r' > N , (A7) 



Q 
f b, 

l Q,N e ,S 

and using the unitary transformation (|A5|) we have, 

M££ S (P> Q)=n<Q, N e , S, p\\f\\Q-l,N e -l,S±^q > N 
= E U &N e s(p,ri)U%_ 1N s± dq,r'i') 

i ■/ 

rxjr'v 

x N ^<Q,N e ,S,ri\\p\\Q-l,N e -hS±^r'i'> N ^ (A8) 
Mj'J e)5 (p, g) = N < Q, N e , S, pp\\Q - 1, N e , S, q > N 

= E u% iN Ap,ri)u»_ 1>Ne Aqyi') 

ri,r'i' 

x N _ x < Q,N e ,S,ri\\tf\\Q-l,N e ,Syi' > JV _ 1 . (A9) 



Evaluating these using the expressions (|Alj-A^) gives an expression relating the matrix 



elements at iteration N to those at iteration N — 1, 



rr' 

+ C£ E KnAPi r2)U»_ x ^ 1}S±l (q, r'2)M££:l s _, (r, r>) 

rr' 

+ C£ E KnAPi rtWl^s^ (q, r'm f Q % N J, s+h (r, O 

rr' 

+ C£ E £$*r.*(P, r^t^L^^i (?, r'4)M^f -; 5 (r, r>) 

rr' 
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and, 



+ Ct 3 E USwfa r2)U%_ hNe _ hS± i (q, r'3)M^-; 5± (r, r') 

rr' 



M b Q% e , s (p, q) = E ^ksfo rl)C/^ liJVej5 (g, rl)M^>, r') 

rr' 

+ E tf&r.s(P, r2)Ug_ w (q, r'2)M b ^ s _ h (r, r') 

rr' 

+ E ^(P. r3)U»_ ljNe>s (q, r'3)M^ 1)5+ , (r, r' 

rr' 

+ E ^ksfe ^Q-i^sfa r'4)M^ 2iS (r, r') 



rr' 



where the coefficients C u , are given in Table |I|. 



APPENDIX B: ERRORS AND SUM RULES 

In this appendix we outline some of the checks carried out to ensure the correctness of 
the numerical renormalization group programs. The eigenvectors and eigenvalues of H can 
be calculated analytically and these can then be used to set up the matrices for N — 1. The 
latter have been compared with those generated by the computer programs for N — 1 and 
found to be identical. We have also checked the recursive evaluation of the matrix elements, 
< r\f^\s >, < r^ls > and < rlcjyjs > by making use of the commutation relations for the 
creation and annihilation operators appearing in the Hamiltonian (HI), 

faft + flU = 1, (Bl) 

Kbl - b% = 1, (B2) 

CNaC ] Na + C ] Na C Na = 1. (B3) 

For any state \Q, N e , S, S z , k >n in the Hilbert space of Hn, the completeness relation 
J2 r \ r >< r \ = 1 together with (|BTHB3"1) yields, 

1 = E I < Qi N e,S,k\\f\\Q -l,N e -l,S -^,k' > | 2 
k> 
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+ ^T2EI <Q,N e ,S,k\\ft\\Q-l,N e -l,S+±V> | 2 
fc' 

+ £| <Q + l,N e + l,S+l-,k'\\p\\Q,N e ,S,k> | 2 , (B4) 
fc' z 

1 = EI <Q+hN e ,S,k'\\bmQ,N e ,S,k> | 2 
fc' 

-531 <Q,iV e ,5,fc||6t||Q-l,iV e ,5,fc'> | 2 , (B5) 
fe' 

1 = X) I <Q,iV e ,^fc||cj v ||Q,iV e -l,S-i fc'> | 2 
fc' 2 

+ tt^EI <g,iv e ,s,fc||c] v ||g,iv e -i, 1 s + ^fc'> | 2 

ZD + Z fc/ Z 

+ £| <Q,iV e + l,5 + ^'||cj v ||Q,iV e ,S , ,fc> | 2 . (B6) 

In the calculations we verified that these relations were satisfied to within rounding errors for 
each state \Q, N e , S,k >iv in for AT = 0, 1, . . . , 4, when all states are retained. This gives a 
reliable test of the formulae in Appendix A and on the routines for the recursive evaluation of 
the matrix elements. We note that once higher energy states start being eliminated, typically 
for A^ > 4, the LHS of the above expressions will be less than unity due to the missing states. 
The sum rule ( |B~6| ) also provides a check on the construction of the Hamiltonian Hn, since 
the latter depends on the matrix elements < Q = 1, N e + 1, S ± |, k'\ \c' N \ \Q = 1, N e , S,k > 



14| . The above tests on all quantities appearing in the iterative procedure for the first few 



iterations virtually eliminates the possibility of errors. 



APPENDIX C: EXACT PROJECTION ONTO THE PHYSICAL SUBSPACE 

In order to effect the constraint of the dynamics to the physical Hilbert subspace it is 
convenient to add the term XQ to the Hamiltonian, where —A is a "chemical potential" 
associated with the auxiliary particle number Q. The operator constraint Q = 1 is imposed 
exactly on the expectation value of any operator O by differentiating with respect to the 
fugacity ( = e~ l3X and then taking the limit A — ► oo 0: 

&ti[Oe-K H+xc »] 

< Q > g = te V^ (g+Aq v (C1) 
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where the trace is taken over the complete, enlarged Hilbert space. In particular, we state the 
following two results which are of use to us in this paper and which follow straightforwardly 
from the above (for details see 0). First, the canonical partition function in the subspace 
Q = 1 is 

Z c = lim tr[Qe-^ +A ( Q - 1 »] 

A^oo 

= \im(e^(Q) G c(X))Z Q=0 , (C2) 

A— >oo 

where the subscripts GC and C denote the grand-canonical and the canonical (Q = 1) 
expectation value, respectively. Second, the canonical Q — 1 expectation value of any 
operator O having a zero expectation value in the Q = subspace is given by, 

(0)c = hm — — — - (C3) 

A^oo (Q) GC {\) 

Thus, we obtain the constrained <i-electron Green function in terms of the grand- 
canonical one (Gd(w,T,\)) as 

G d {u) = hm — - — — (C4) 

A^oo {Q) GC {\) 

In the enlarged Hilbert space (Q = 0, 1, 2, ...) Gd(u, T, A) may be expressed in terms of the 
pseudo-fermion and slave boson Green functions using Wick's theorem. It then follows from 
Eq. flC4j) that the operator constraint Q = 1 is imposed on the auxiliary Green's functions 
by simply taking the limit A — > oo of the respective unconstrained functions. Clearly, by 
this procedure all excitation energies of pseudo-fermions and slave bosons are shifted to oo. 
It is therefore convenient to re-define the auxiliary particle frequency scale as u — > u + A 
before taking the limit A — ► oo. Note that this does not affect the energy scale of physical 
quantities (like the local d electron Green's function), which is the difference between the 
the pseudo-fermion and the slave-boson energy. 
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APPENDIX D: NCA CALCULATION OF THE AUXILIARY SPECTRAL 

FUNCTIONS Ap B 

In order to enter the asymptotic power law regime of the auxiliary spectral functions and 
to compare with the T=0 results of the NRG, the NCA must be evaluated for temperatures 
several orders of magnitude below T , the low temperature scale of the model. The equations 
are solved numerically by iteration. In this appendix the two main procedures are described 
to make the diagrammatic auxiliary particle technique suitable for the lowest temperatures. 

The grand-canonical expectation value of the auxiliary particle number appearing 
in Eq. ( |CJ4| ) is given in terms of the (unprojected) auxiliary particle spectral functions 
A+ b (uj,T, A) by 

(Q)gc(X) = I duj[f(uj)J2AU^^,X) + b(uj)AUoo,T 1 X)} (Dl) 

where f(ui), b(u>) denote the Fermi and Bose functions. Substituting this into the expression 
( |U2[ ) for the canonical partition function we obtain after carrying out the transformation 
uj — > uj + A, and taking the limit A — > oo 

e-^CO = = lim e^(Q) GC (X) 

= / du;e-^[£A%(u;,T)+A+(u;,T)}, (D2) 

where A^ b (uj,T) are now the projected spectral functions as defined in Eq. @, and by 
definition F imp = — ^\ii(Zc/Zq = q) is the impurity contribution to the Free energy. 

The numerical evaluation of expectation values like (Q)gc(X — > oo) (Eq. (P2|) ) or 
E CCT (a;,T, A — > oo) (Eq. (p2])) is non-trivial (1) because at T = the auxiliary spectral 
functions A~f b (u,T) are divergent at the threshold frequency Eq, where the exact position 
of E is a priori not known, and (2) because the Boltzmann factors e _/3w diverge strongly 
for uj < 0. Therefore, we apply the following transformations: 

(1) Before performing the projection uj uj + X, A ^ oo we re-define the frequency 
scale of all auxiliary particle functions A^ b according to uj — > uj + A , where A is a finite 
parameter. In each iteration Ao is then determined such that 
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J<\ue-^A%{u) + At{u)] = l (D3) 

where A^ b {uj) = Hindoo A~f b (uj + A + A, T, A) is now an auxiliary spectral function with the 
new reference energy. It is seen from Eq. ( |D2| ) that A (T) = F imp (T) = Fq = %(T) —Fq =0 (T), 
i.e. Ao is the chemical potential for the auxiliary particle number Q, or equivalently the 
impurity contribution to the Free energy. The difference of the Free energies becomes equal 
to the threshold energy E = E%l x - Egl at T = 0, so the energy scale of the shifted 
spectral functions defined above coincides exactly with that of the NRG spectral functions 
defined in (||). More importantly, however, the above way of determining a "threshold" is 
less ad hoc than, for example, defining it by a maximum in some function appearing in the 
NCA equations. It is also seen from Eq. ([D3|) that this procedure defines the frequency scale 
of the auxiliary particles such that the T = threshold divergence of the spectral functions 
is at the fixed frequency u = 0. This substantially increases the precision as well as the speed 
of numerical evaluations. Eq. ( p4j ) for the projected d electron Green's function becomes 

G d (u) = \ime^G d (uj,T,X). (D4) 

(2) The divergence of the Boltzmann factors implies that the self-consistent solutions 
for Af b (u) vanish exponentially ~ e^ 3 for negative frequencies. It is convenient, not to 
formulate the self-consistent equations in terms of A^ b like in earlier evaluations ||, but to 
define new functions Af^(u) and ImEf^u) such that 

4*>H = /(-") A f^) ( D5 ) 

Im£/ i6 (a;) = f{—uS) ImS/ i6 (o;). (D6) 

After fixing the chemical potential Ao and performing the projection onto the physical sub- 
space, the canonical partition function (Eq. ( |C2|) ) behaves as liniA^oo e™ x ~ Xo > Zc(T) = 1, 
and it follows immediately from the definition of Aj b that 

AJ^lj) = f(u) A f , b {u). (D7) 
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In this way all exponential divergencies are absorbed by one single function for each particle 
species. The NCA equations in terms of these functions are free of divergencies of the 
statistical factors and read 

ImE /ff (u, - ,0, T) = V 2 Y / (l ~ /(£fc 1 ))(l : / ^~ £fc)) i 6 (a; - e k ) (D8) 
ImE^ - i0, T) = V 2 Y, f{€k){ )~{ { ^ ek)) A fa ^ + e k ), (D9) 

(g)(A ,A^oo)= /d w /HEi /(7 H + 4H] = i (dio) 

ImGW(^ - ^0, T) — [ de[f(e + u)f(-e) + /(-e - u)f(e)]A fa (e + co)A b (e), (Dll) 



where the real parts of the self-energies E/, E&, E c are determined from a Kramers-Kronig 
relation, and the auxiliary spectral functions are the imaginary parts of the Green's functions, 



A%(u; 



— Im 



,A+(u) = -Im 



{uj + Xo-iO-^biuj-iO))' 1 



u+\o-iO-e d -E fcr (u)-iO)) 1 
The above method allows to solve the NCA equations effectively for temperatures down 
to typically T = 10 _4 To. It may be shown that the same procedure can also be applied to 
self-consistently compute vertex corrections |E3| beyond the NCA. 
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FIGURES 

FIG. 1. (a) The lowest order contribution to the generating functional <F and, (b), the renor- 
malized vertex part A(ico u , iu> u —iio), entering the expression ( p6| ) for the d-electron Green's function 
(with frequency conservation we can omit iu n in A). The solid lines are for band electrons, the 
dashed lines are for pseudo-fermions and the wiggly lines are for slave-bosons. 

FIG. 2. The T = NRG pseudo-fermion (solid lines) and slave-boson (dashed lines) 
spectral functions in the Kondo case e^/A = —4, To/A = 1.87 x 10~ 2 , = 0.874. The + signs 
are for the spectral function above the threshold, Eq, and the circles are for the spectral function 
below the threshold. The arrow indicates the position of \ed\- 

FIG. 3. The T = NRG spectral functions A^ b in the mixed valent regime e^/A = 0, 
= 0.314 with notation as in Fig. |2[ 

FIG. 4. The T = NRG spectral functions A^ b in the empty orbital regime e^/A = +2, 
rid = 0.172, with notation as in Fig. B 

FIG. 5. The exponents aj (o), a& (o) deduced from the asymptotic power law behaviour of 
the auxiliary spectral functions as calculated within the NRG for different values of the occupation 
rid- The solid lines are the functions — n\/2 and 1 — tVj/2. 

FIG. 6. The temperature dependence of the occupation number n^(T) for different (the 
curves are labeled by e^/A). The high temperature limit of 2/3 (indicated by an arrow) is recovered 
in all cases. 

FIG. 7. Comparison of the NRG (o) and NCA (o) auxiliary spectral functions A^ b (a-b), Aj b 
(c-d) in the Kondo case e^/A = —4. The arrow indicates the position of |e^|. The NRG results 
are for T = and the NCA results are for T = 1.0 x 10~ 6 D. The diver gence of the NCA spectral 
functions for to — > is cut off below the finite temperature T << Tq. 
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FIG. 8. Comparison of the NRG (o) and NCA (o) auxiliary spectral functions A^ b (a— b), A^ b 
(c-d) in the mixed valent case e^/A = 0. The NRG results are for T = and the NCA results are 
for T = 1.0 x 10~ Q D. 

FIG. 9. Comparison of the NRG (o) and NCA (o) auxiliary spectral functions A^ b (a— b), Aj b 
(c-d) in the empty orbital case e^/A = +2. The NRG results are for T = and the NCA results 
are for T = 1.0 x 10" 6 D. 

FIG. 10. The impurity spectral function in the Kondo regime = — 4A, To/A = 1.87 x 10 -2 , 
for (a) high energies and, the low energy region (b) —5 < lo/Tq < 10, and (c) —2 < lu/Tq < +2. 
The dashed curve is the NRG result p% RG ', the dot-dashed curve is the NRG result without the 
vertex part in (|33[), and the solid curve is the NCA result. The NRG results are for T = and the 
NCA results are for T = 1.0 x 10- 6 D. 

FIG. 11. The impurity spectral function in the mixed valent regime e^/A = 0, with notation 
as in Fig. [TO. The NRG results are for T = and the NCA results are for T = 1.0 X W~ 6 D. 



FIG. 12. The impurity spectral function in the empty orbital regime e^/A = +2, with notation 
as in Fig. M. The NRG results are for T = and the NCA results are for T = 1.0 X 10" 6 D. 
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TABLES 

TABLE I. The threshold exponents aty, for the auxiliary spectral functions 
A^ b ~ \uj — Eo\~ a f' h . The quantities a'^ b are n^ — n d /2 and 1 — n^/2, respectively where rid = n^ RG 
is the impurity occupation calculated from the NRG partition function (quantities shown to 3 signif- 
icant figures). The NCA results for the impurity occupation, n^ CA , are also shown. The impurity 
spectral density at the Fermi level calculated within the NCA, p^ CA (eF), and NRG, p^ RG (ep), 
are tabulated together with the % relative error in the Friedel sum rule, Pd{^F) = sin 2 (7rn(i/2)/7rA. 
The exponent ctf is difficult to estimate close to = 1 due to the small Kondo scale and the slow 
asymptotic behaviour of the pseudo-fermion spectral function. The low energy scale To is Tx given 
by ( |27| ) in the Kondo regime, A in the mixed valent regime and in the empty orbital regime. 
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0.498 


0.563 


0.564 


-5 
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TABLE II. Coefficients C% 



a' c% cz 

IV IV 

11 1 1 



99 / 25(25+2) -, 

*' (25+l)(25+l) 1 



25(25+2) 
(25+l)(25+l) 



33 1 
44 1 1 

23 25TT 

32 _ 25TT 
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